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I.  INTRODUCTION 


This  project  explores  the  uses  of  time  domain  capable  vector 
network  analyzers  as  a  diagnostic  in  helix  traveling-wave  tubes.  The 
automatic  network  analyzer  (ANA)  used  throughout  the  experiments  was  the 
Hewlett  Packard  8510  ANA  with  the  time  domain  option. 

In  Section  II  of  this  paper,  the  experimental  set-up  used  for 
these  experiments  is  discussed. 

In  Section  III  of  the  paper,  the  uses  of  the  time  domain  are 
investigated.  The  first  section  gives  some  general  theory  of  time 
domain  analysis  and  some  of  the  items  of  concern  when  using  a  network 
analyzer  with  time  domain  capability.  The  preliminary  investigations 
involved  looking  at  the  typical  response  of  helix  circuits,  both  for 
special  cold  test  circuits  and  actual  helix  circuits.  Also  investigated 
was  the  ability  of  the  ANA  to  resolve  a  bead  perturber  as  it  moved  along 
the  circuit.  Next,  a  method  of  determining  the  phase  velocity  versus 
frequency  was  devised  and  tested  on  ideal  circuits,  an  airline  circuit, 
and  finally  on  actual  tube  circuits.  The  final  use  of  the  time  domain 
that  was  investigated  was  the  more  traditional  use  of  locating  faults  or 
discontinuities  as  a  function  of  distance  in  the  circuits. 

In  Section  IV  of  the  paper,  the  complex  vector  nature  of  the  data 
obtained  from  the  ANA  is  exploited  to  enable  one  to  obtain  phase  veloc¬ 
ity  as  a  function  of  frequency  without  going  into  the  time  domain.  The 
general  theory  behind  the  method  is  discussed,  along  with  the  model 
used  to  calculate  the  phase  velocity  and  the  method  of  taking  data.  The 
results  are  then  compared  with  previously  obtained  data  to  illustrate 
the  accuracy  of  the  method. 
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II.  EXPERIMENTAL  SET-UP 


The  measurements  in  this  paper  were  all  made  using  a  metallic  bead 
as  a  perturbing  object,  as  has  been  suggested  in  the  literature  dealing 
with  perturbation  techniques. 1  The  bead  was  suspended  on  a  thin  nylon 
line  and  drawn  down  the  axis  of  the  circuit.  In  these  types  of  measure¬ 
ments,  it  is  imperative  to  insure  that  the  perturbing  object  is  centered 
in  the  circuit  being  measured  in  order  to  get  accurate  results. 

In  order  to  get  some  control  over'  the  position  of  the  bead  in  the 
circuits,  a  fixture  was  made  to  hold  the  circuits  while  measurements 
were  being  made.  To  facilitate  centering  of  the  bead,  the  circuit  was 
placed  in  a  cradle  and  secured  in  place  with  tie  downs.  This  cradle  was 
designed  to  have  some  movement  in  the  horizontal  direction  so  as  to  be 
able  to  get  the  string  straight  down  the  circuit.  Vertical  adjustabil¬ 
ity  was  incorporated  into  the  connection  to  the  slotted  line  carriage 
and  also  to  the  pulley  at  the  other  end.  The  line  is  attached  to  the 
carriage,  threaded  through  the  circuit,  and  then  hung  over  the  pulley 
with  a  weight  to  provide  tension  in  the  line.  With  adjustability  in 
both  directions,  it  is  then  easy  to  get*  the  bead  to  travel  along  the 
axis.  The  slotted  line  carriage  has  a  scale  on  it  that  is  accurate  down 
to  one  hundred  micrometers.  For  increased  accuracy  in  positioning,  a 
high  resolution  stepper  motor  would  be  the  best  solution,  and  should 
result  in  much  better  repeatability  than  moving  the  carriage  by  hand.  A 
stepper  motor  could  also  be  programmed  to  make  the  movement  of  the  bead 


automatic. 


The  ANA  was  used  both  manually  and  controlled  by  an  HP  series  200 
computer.  The  computer  could  be  used  to  control  a  stepper  motor  to  make 
the  measurements  semiautomatic,  although  some  user  input  would  still  be 
needed  for  the  time  domain  methods. 


III.  TIME  DOMAIN  MEASUREMENTS 

A.  General  Theory  and  Capabilities 

Before  going  into  the  results  of  the  investigation  of  the  uses  of 
the  TD  capabilities  for  TWT's,  an  overview  of  the  general  theory  and  the 
specific  abilities  of  the  HP8510  vector  network  analyzer  will  be  pre¬ 
sented.  For  a  more  in-depth  description  of  time  domain  theory,  the 
reader  is  referred  to  the  literature2  and  to  the  network  analyzer  manu¬ 
als  for’  a  more  detailed  description  of  the  ANA’s  abilities. 

Basically,  the  problem  is  to  take  complex  frequency  domain  data 
and  transform  it  to  the  time  domain  using  some  form  of  a  Fourier  trans¬ 
form.  The  relationship  between  the  frequency  domain  response  and  the 
time  domain  response  is  described  by  the  Fourier  transform: 

Frequency  domain  Time  domain 
H(f) - ->•  h(T) 

Coupling  this  information  with  the  frequency  response  of  the  device  as  a 
function  of  frequency,  one  can  obtain  information  on  the  presence,  and 
sometimes  type,  of  mismatches  in  the  device  as  a  function  of  distance. 
This  can  be  a  powerful  aid  in  determining  the  location  of  faults  and 


also  in  seeing  where  the  various  discontinuities  are  located  in  the 
circuit.  The  important  thing  to  remember  is  that  the  frequency  domain 
measurement  is  a  composite  response  of  all  the  discontinuities  in  the 
device,  while  the  time  domain  measurement  shows  the  effect  of  each 
individual  discontinuity  as  a  function  of  time  and  distance. 
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There  are  two  modes  of  operation  in  the  time  domain  when  using  the 
HP8510  ANA.  The  time  domain  bandpass  mode  is  useful  for  measuring  band- 
limited  devices,  while  the  time  domain  lowpass  mode  mathematically 
simulates  the  traditional  type  of  TDR  measurement  and  gives  the  user 
information  on  the  type  of  discontinuity  present,  as  well  as  its  physi¬ 
cal  location  in  the  device  under  test. 

The  bandpass  mode  has  the  advantage  that  there  is  no  restriction 
on  the  frequency  span  of  the  measurement,  although  resolution  and  other 
factors  can  put  a  practical  limit  on  the  bandwidth  of  the  measurement. 
The  lowpass  mode,  on  the  other  hand,  is  limited  in  its  use,  since  the 
nature  of  the  transformation  requires  that  frequency  data  be  taken  all 
the  way  down  to  zero  frequency,  and  thus  restricts  its  use  with  bandpass 
devices.  This  mode  gives  the  user  some  information  that  can  be  used  to 
determine  the  type  of  discontinuity  that  is  present,  and  also  has  the 
best  resolution  of  the  two  modes.  Some  concepts  of  interest  to  any  user 
of  time  domain  measurement  systems  include  the  effects  of  masking, 
windowing  and,  most  importantly,  resolution.  Masking  is  a  physical 
phenomenon  in  which  the  response  of  one  discontinuity  affects  the 
response  of  all  the  subsequent  discontinuities  in  the  circuit.  This  is 
a  result  of  the  fact  that  the  power  lost  at  each  discontinuity  will 
never  reach  the  mismatches  further  on  in  the  circuit.  Windowing  is  a 
technique  used  to  try  to  limit  the  adverse  effects  of  the  abrupt  transi¬ 
tion  in  the  data  at  the  start  and  stop  frequencies.  This  transition 
introduces  undesirable  effects,  such  as  overshoot  and  ringing  in  the 
time  domain  response.  When  this  type  of  data  is  transformed  to  the  time 
domain,  it  has  a  sin(kt)/kt  shape,  where  k  is  a  constant  and  t  is 
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time.  The  resultant  finite  impulse  width  limits  the  measurement’s 
resolution  and  sidelobes  are  introduced  that  further  limit  the  dynamic 
range  of  the  measurement. 

Windowing  essentially  filters  the  data  in  the  frequency  domain  to 
smooth  out  the  abrupt  transitions  and,  thus,  can  reduce  the  aforemen¬ 
tioned  deleterious  effects.  By  varying  the  shape  of  the  window,  one  can 
enhance  the  data,  but  a  tradeoff  is  made.  For  example,  sidelobes  can  be 
reduced,  but  a  price  is  paid  in  increased  impulse  width. 

The  most  important  considerations  are  range  and  resolution.  Range 
is  the  amount  of  distance  one  can  see  into  the  device  before  encounter¬ 
ing  the  aliasing  effects  that  occur  for  this  type  of  discrete  transfor¬ 
mation.  In  other  words,  the  nature  of  the  transformation  causes  the 
response  to  be  repeated  at  regular  intervals.  Range  is  primarily  depen¬ 
dent  on  the  frequency  spacing  of  the  data  points.  There  are  two  kinds 
of  resolution  to  be  concerned  with:  range  resolution  and  response 
resolution.  Range  resolution  is  the  ability  to  locate  a  mismatch  in 
time,  or  how  accurately  its  position  can  be  pinpointed.  Response  reso¬ 
lution  is  how  far  apart  two  separate  discontinuities  have  to  be  in  order 
for  one  to  be  able  to  see  them  distinctly.  The  response  resolution  is 
of  most  importance  and  is  a  function  of  the  window  shape  and  also 
inversely  proportional  to  the  frequency  span  of  the  measurement. 

The  feature  of  the  transformation  that  is  of  greatest  interest  is 
the  ability  to  filter  the  data  in  the  time  domain  to  isolate  an  individ¬ 
ual  response,  and  then  transform  back  to  the  frequency  domain  to  view 
the  frequency  response  of  only  that  response.  This  is  a  powerful  tool, 
since  theoretically  one  could  look  inside  a  device  and  observe  the 


frequency  response  of  each  individual  response  instead  of  the  composite 


response. 

One  important  aspect  of  looking  at  the  time  domain  response  of  a 
circuit  is  that  the  response  is  necessarily  an  average  over  the  entire 
bandwidth  in  the  frequency  domain.  For  example,  any  of  the  responses 
that  make  up  the  composite  response  for  a  device  in  the  frequency  domain 
can  be  singled  out  in  the  time  domain,  but  any  data  taken  in  this  mode 
would  be  an  average  over  the  frequency  span.  This  is  an  important 
limitation,  since  one  usually  wishes  to  know  the  phase  velocity  as  a 
function  of  frequency  in  the  band  of  interest. 

B.  Preliminary  Investigation 

1 .  Response  of  Helix  Circuits 
a.  Cold  Test  fixture  Versus  Actual  Circuits 

The  experiments  using  the  time  domain  capabilities  of  the  HP8510 
ANA  were  done  on  both  actual  helix  circuits  and  special  cold  test  fix¬ 
tures.  The  cold  test  fixtures  are  helix  circuits  that  are  specially 
assembled  for  doing  tests,  such  as  phase  velocity.  Typically,  the 
support  rods  in  such  circuits  have  no  loss  deposited  on  them,  and  the 
helix  itself  might  be  scaled  up  in  diameter  in  order  to  facilitate 
tests,  such  as  pulling  a  bead  through  the  circuit.  In  addition,  the 
helix  is  usually  of  uniform  pitch.  The  helix  circuits  tested,  on  the 
other  hand,  are  circuits  that  are  taken  from  the  production  line  and 
tested  fn  situ.  The  rods  have  a  loss  pattern  deposited  on  them  and  the 
helix  might  have  a  pitch  taper  at  some  point  in  the  circuit.  Most  of 
the  work  was  done  on  actual  circuits;  however,  the  work  done  on  cold 
test  circuits  is  also  discussed. 


-  7  - 


u.  lypicai  nesuuuae 


The  first  thing  done  was  to  get  some  familiarity  with  what  the 
time  domain  response  of  helix  circuits  looked  like.  A  typical  response 
is  shown  in  Fig.  1,  where  the  two  large  discontinuities  are  composite 
discontinuities  that  are  made  up  of  the  connector  and  pin-to-helix 
connection.  As  can  be  seen  in  the  figure,  the  area  between  these  two 
discontinuities  is  essentially  flat  indicating  no  substantial  mismatches 
along  the  helix  length.  When  a  bead  supported  on  a  length  of  nylon 
fishing  line  was  inserted  into  the  circuit,  a  typical  response  is  shown 
in  Fig.  2.  This  figure  gives  some  idea  of  the  relative  magnitude  of  the 
bead  response  as  compared  to  the  responses  of  the  input  and  output.  The 
bead  was  approximately  three-fourths  of  the  helix  diameter  and  one  pitch 
distance  wide,  and  was  made  by  either  putting  a  small  amount  of  silver 
loaded  epoxy  on  the  nylon  line  or  by  gluing  a  3mall  piece  of  metal 
tubing  to  the  line. 


2.  Bead  Tracking 


The  first  test  run  was  to  see  if  the  movement  of  the  bead  down  the 
circuit  could  be  tracked  accurately.  The  bead  was  moved  a  known  dis¬ 
tance,  and  then  the  distance  from  the  reference  plane  to  the  peak  of  the 
bead  response  was  measured  on  the  ANA.  As  is  typical  with  time  domain 
reflection  measurements,  the  distance  is  the  round  trip  distance  and 
must  also  be  corrected  for  the  circuit  phase  velocity.  The  relative 
phase  velocity  (vp/c)  for  this  cold  test  circuit  was  0.205.  Table  1 
gives  the  results  of  the  tests. 
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Fig.  1.  A  typical  response  of  a  helix  cold  test  circuit. 

The  first  response  is  due  to  the  input  disconti¬ 
nuity,  and  the  second  to  the  output  discontinuity. 
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Fig.  2.  The  response  of  a  helix  cold  test  fixture  with  no 
bead,  and  with  the  bead  at  two  different  positions 
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Table  1.  Comparison  of  the  physical  distance  moved  by  the  bead 
(z  column)  with  the  value  calculated  from  the  peak  of 
the  bead’s  response  in  the  time  domain  (delta  z  column) 
v/c  was  assumed  to  be  0.205  to  make  the  calculations. 


z(cm) 


Delta  z(cm) 


0.6038 


0.5531 


0. 502*1 


0.5623 

0.5670 

0.5716 

0.3*157 

0.1179*1 

0.11563 


Percent 

Difference 


20.76 


10.62 


12. H6 


13.H0 

1*4.31 

-30.86 

-4.126 

-8.736 


These  results  were  not  very  encouraging,  so  another  method  was 
used.  It  was  possible  that  the  true  position  or  the  bead  was  not  being 
accurately  found.  One  way  to  try  and  pinpoint  the  position  is  to  memo¬ 
rize  the  trace  in  the  time  domain  response  of  the  circuit  without  the 
bead  present,  and  then  subtract  this  response  from  the  response  when  the 
bead  is  present.  This  way  the  ANA  should  show  more  accurately  the 
difference  that  the  bead  is  making  to  the  response.  However,  the  accu¬ 
racy  in  the  ANA's  ability  to  track  the  bead  was  still  poor,  as  is  seen 


Jf 
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in  Table  2. 


Table  2.  Comparison  of  the  actual  bead  location  to 
peak-of-bead  response  in  the  time  domain. 


z(cm) 

Marker  Value 
(cm) 

0.0 

1.0052 

1.0 

1.7245 

2.0 

2.7943 

3.0 

3.9748 

4.0 

5.3120 

Delta  z 
(cm) 


0.7193 

1.0698 

1.1805 

1.3373 


Percent 

Difference 


-28.1 

6.98 

18.05 

33.73 


'.  Determination  of  vp  Versus  f  Using  Gatin 


In  this  section,  the  use  of  gating  as  a  tool  to  measure  the  phase 
velocity  is  investigated.  The  method  of  measurement  is  to  insert  the 
bead  into  the  circuit  and  then  transform  it  to  the  time  domain.  The 
response  of  the  bead  alone  can  be  gated  upon  and  then  transformed  back 
to  the  frequency  domain.  This  will  give  the  complex  reflection  coeffi¬ 
cient  as  a  function  of  frequency  of  the  bead  referenced  to  the  reference 
plane  of  the  network  analyzer.  The  value  of  the  bead’s  phase  angle  is 
noted,  and  the  bead  is  then  shifted  a  known  distance  to  a  new  position 
where  the  process  is  repeated.  The  resultant  phase  shift  should  corre¬ 
late  with  the  circuits  phase  velocity  at  each  frequency. 

This  method  depends  on  the  ability  of  the  network  analyzer  to 
precisely  locate  the  bead  within  the  circuit.  As  has  already  been 
discussed,  the  ANA  is  not  accurate  in  this  regard.  Also  of  concern  here 
is  the  quality  of  the  filtering  done  in  the  gating  process.  Since  the 
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bead's  response  is  typically  small  compared  with  the  response  of  the 
input  discontinuity ,  if  the  filter  is  not  narrow  enough,  some  of  the 
input  discontinuities  response  will  be  included  and  will  affect  the 
accuracy  of  the  data  in  the  frequency  domain.  This  is  especially  a 
problem  at  the  band  edges,  although  the  data  towards  midband  is  usually 
unaffected  unless  the  bead  is  very  close  to  the  input  discontinuity. 

As  a  check  to  see  if  this  method  has  any  merit,  some  tests  were 
done  on  ideal  circuits  using  the  circuit  modeling  program  (CMP)  that  is 
included  with  the  HP8510  ANA  time  domain  option,  and  also  on  some 
slotted  airline. 


1.  Verification  of  Method 


a.  Circuit  Modeling  Program  (Ideal  Case) 

The  circuit  modeling  program  is  a  program  which  can  be  used  to 
compute  the  s  parameters  of  ideal  circuits  made  up  of  transmission 
lines,  resistors,  capacitors,  and  inductors.  This  information  is  then 
fed  into  the  HP8510  and  treated  as  normal  data.  It  can  then  be  operated 
on  and  displayed  in  any  manner  that  the  ANA  allows.  In  this  case, 
typical  responses  seen  in  the  helix  circuits  were  modeled  using  trans¬ 
mission  lines,  and  then  manipulated  in  order  to  get  the  frequency 
response  of  the  modeled  bead  alone.  This  allows  verification  of  the 
validity  of  the  gating  method  for  determining  phase  velocity.  Also,  it 
was  used  to  observe  the  effect  of  composite  discontinuities,  such  as 
those  which  appear  in  actual  circuits  at  the  input  as  well  as  what 
happens  when  the  phase  velocity  changes  along  the  circuit. 
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The  first  verification  test  was  a  simple  cascade  of  a  50  ohm 
transmission  line  into  a  75  ohm  transmission  line  and  back  into  a  50  ohm 


line.  9.99  GHz  was  arbitrarily  chosen  as  the  frequency  at  which  to 
verify  the  method.  This  results  in  two  mismatches  in  the  time  domain. 
What  is  of  interest  here  is  to  see  if  changing  the  length  of  the  75  ohm 
line  results  in  the  expected  behavior  of  the  angle  of  the  complex 
reflection  coefficient  of  the  second  mismatch.  Since  the  phase  velocity 
is  the  speed  of  light,  we  have  at  9.99  GHz, 

0  -  2*ir*9.99  GHz/2. 99e9  »  209.23  rad/meter  *  H9.88°/cm 

The  results  of  measuring  the  total  reflection,  gated  first  mismatch,  and 
gated  second  mismatch  angles  are  shown  in  Table  3.  The  expected  shift 
in  phase  20^  for  a  0.2  cm  shift  in  the  length  of  the  75  0  line  is  ^7 - 952 
degrees,  and  119.88  degrees  for  a  shift  of  0.5  cm,  using  the  0  calcu¬ 
lated  above.  From  Table  3,  it  can  be  seen  that  the  shift  ir.  phase  was 
107.25  -  59.25  -  47.992  for  0.2  cm  and  179.22  -  107.25  -  119.96  for  0.5 
cm.  These  results  are  well  within  a  percent  of  the  expected  values! 

Table  3.  The  angle  of  the  reflection  coefficient  for  the 
total  circuit,  the  first  mismatch  alone,  and  the 
second  mismatch  alone  for  an  ideal  circuit  con¬ 
sisting  of  a  50  to  75  to  50  ohm  transmission  line. 


L(cm) 

Total 

Reflection 

Angle 

First 

Reflection 

Angle 

Second 

Reflection 

Angle 

10.0 

-27.701 

0.0013733 

-59.258 

10.2 

-51. *11 

0.003733 

-107.25 

10.5 

-89.58 

0.004H99 

-179.22 
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As  stated  above,  it  is  also  of  interest  to  see  if  this  method  is 
accurate  when  the  types  of  composite  discontinuities  seen  in  helix 
circuits  are  present.  The  response  of  a  bead  in  a  helix  circuit  was 
modeled  by  placing  a  composite  mismatch  in  front  of  the  response  of  the 
bead.  It  was  attempted  to  keep  the  relative  magnitudes  close  to  what  is 
actually  seen  in  helix  circuits.  The  time  domain  response  of  the  model 
is  shown  in  Fig.  3a,  and  the  equivalent  circuit  shown  in  Fig.  3b.  The 
length  of  the  76  ohm  TL  (simulating  the  shifting  of  the  bead's  response 
as  it  is  moved  down  the  circuit)  was  increased  0.5  cm,  and  the  phase  of 
the  gated  reflection  at  9.99  GHz  changed  by  120  degrees  as  expected. 
This  indicates  that  the  composite  response  does  not  adversely  affect  the 
gating  method  in  the  ideal  case.  The  last  question  to  be  addressed  is 
the  effect  of  changing  phase  velocity  on  the  method.  In  an  actual 
circuit,  one  would  expect  the  region  from  the  input  window  to  the  pin¬ 
to-helix  connection  to  have  a  phase  velocity  roughly  corresponding  to 
coaxial  line.  The  region  after  the  input  discontinuity  will  have  a 
phase  velocity  characteristic  of  the  helix  itself.  For  this  reason, 
another  test  case  was  run  with  an  equivalent  circuit  as  shown  in  Fig. 
*4.  The  length  of  70.5  ohm  TL  was  set  to  a  phase  velocity  corresponding 
to  the  speed  of  light,  while  the  82  ohm  line  was  given  a  relative  phase 
velocity  (v/c)  of  0.178,  which  is  a  reasonable  value  for  a  helix  cir¬ 
cuit.  The  beta  of  the  82  ohm  line  is  now  67^.16  decrees/cm  at  9.99  GHz, 
and  thus  one  expects  a  shift  of  67^.16  degrees  (.26  L)  for  a  shift  in 
distance  of  0.5  cm.  The  phase  shift  obtained  by  gating  was  673.9^, 
which  is  within  a  percent  of  the  expected  value  and  shows  that  theoreti¬ 
cally  the  changing  phase  velocity  should  not  affect  the  measurement. 
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Fig.  3a.  The  ideal  response  of  the  circuit  shown  in  Fig.  3b. 

This  response  is  intended  to  model  the  type  of  response 
seen  in  helix  circuits.  The  first  response  is  a  com¬ 
posite  response  due  to  the  input  window/helix-to-pin 
connection,  and  the  last  response  is  due  to  the  bead. 


b .  Airline 

Since  the  method  seems  to  work  properly  in  the  ideal  case,  it  was 
tried  on  a  slotted  airline  fixture.  This  seemed  to  be  the  next  step  up 
in  complexity.  Two  calibrations  were  used:  Calibration  1  was  from  0-22 
GHz  to  allow  for  use  of  the  lowpass  mode,  and  Calibration  2  was  from 
7-17  GHz  and  was  used  in  the  bandpass  mode  of  the  time  domain.  The  gat¬ 
ing  method  as  previously  described  was  used,  and  the  results  are  shown 
in  Tables  4a  and  4b  for  both  calibrations.  Since  airline  essentially 


Table  4a.  The  phase  shift  resulting  from  several  con¬ 
secutive  0.5  cm  shifts  in  the  bead  position 
at  10  GHz.  The  expected  shift  is  120  degrees. 


Delta  z  (cm) 

Cal  1: 

Delta  $ 
(Degrees) 

Cal  2: 

Delta  <f> 
(Degrees) 

0.5 

122.63 

121.99 

0.5 

124.14 

121.27 

0.5 

123.84 

123.62 

0.5 

126.38 

124.49 

Table  4b.  The  Phase  shift  resulting  from  two  consec¬ 
utive  0.5  cm  shifts  in  the  bead  location  at 
12  GHz.  The  expected  shift  is  144  degrees. 


Cal  1: 

Cal  2: 

Delta  $ 

Delta  <(> 

Delta  z  (cm) 

(Degrees) 

(Degrees) 

0.5 

141.19 

138.31 

0.5 

158.82 

162.68 

m 

% 
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has  a  normalized  phase  velocity  of  1.0,  the  expected  phase  changes  for  a 
0.5  cm  shift  at  10  GHz  was  again  120  degrees,  while  at  12  GHz  the 
expected  phase  shift  was  144  degrees.  As  can  be  seen  in  Table  4a , 
agreement  is  within  5  percent  at  10  GHz  for  both  calibrations ,  while 
from  Table  4b  the  data  at  12  GHz  were  not  as  good. 

This  simple  test  was  deemed  encouraging  enough  to  have  some  faith 
in  the  method.  The  test  also  seemed  to  indicate  that  lowpass  mode  and 
bandpass  mode  give  roughly  the  same  results  for  this  type  of  measure¬ 
ment.  The  next  step  was  to  try  the  method  on  some  actual  helix  circuits. 

2.  Testing  of  Actual  Circuits 

The  first  circuit  tested  was  a  7-17  GHz  circuit  with  a  mean  diame¬ 
ter  of  0.0723  inches  and  a  pitch  of  0.045  inches.  The  results  of  test¬ 
ing  this  circuit  were  not  very  encouraging.  The  response  of  the  circuit 
when  the  bead  was  gated  on  did  not  give  any  repeatable  results.  The 
phase  angles  did  not  correlate  with  expected  phase  shifts,  and  also  did 
not  even  shift  the  same  in  angle  for  the  same  shift  in  distance.  The 
next  circuit  tested  was  at  a  lower  frequency  and,  as  a  result,  also  had 
a  larger  diameter,  making  it  easier  to  keep  the  bead  centered  as  it  was 
drawn  down  the  axis.  The  band  for  which  the  circuit  was  designed  was  2- 
8  GHz,  and  the  mean  diameter  and  pitch  were  0.151  inches  and  0.088 
inches,  respectively.  The  bead  used  in  this  case  had  a  diameter  of 
0.110  inches  and  a  width  of  0.054  inches.  The  results  of  the  experiment 
are  3howr.  plotted  in  Fig.  5  against  data  obtained  for  a  cold  test  fix¬ 
ture  made  with  the  same  circuit  and  tested  using  the  resonance  method 
discussed  in  Appendix  C.  The  data  are  also  compared  numerically  in 
Table  5. 
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es 

onance  method  on  a  2-8  GH 

z  c 

ka  cot(ii/) 

Gating 

(v/c)cot(iJ/) 

Resonance 
( v/c)cot(  l|>) 

Percent 

Difference 

0.4334 

1.1205 

0.9467 

18.36 

0.5634 

1.0201 

0.9402 

8.50 

0.6935 

0.9647 

0.9348 

3.20 

0.8235 

0.9327 

0.9274 

0.57 

0.9535 

0.9152 

0.9213 

-0.66 

1.0836 

0.9080 

0.9238 

-1.71 

1.2136 

0.9074 

0.9130 

-0.61 

1.3436 

0.9058 

0.9087 

-0.32 

1.4737 

0.9019 

0.9058 

-0.43 

1.6037 

0.9004 

0.8951 

0.59 

1.7337 

0.8392 

0.8973 

-6.48 

From  the  data  presented,  it  can  be  seen  that  agreement  is  very 
good  in  the  middle  of  the  band,  but  is  not  good  at  the  band  edges.  This 
type  of  degradation  in  the  accuracy  of  the  response  at  the  band  edges  is 
related  to  the  fact  that  the  gating  is  not  perfect.  If  the  gate  is  not 
centered  accurately  or  the  span  is  set  too  wide,  the  frequency  domain 
response  will  not  be  the  response  of  the  gated  mismatch  alone,  but  will 
also  be  affected  by  the  other  mismatches  in  the  circuit.  In  any  case, 
some  of  the  effects  of  the  rest  of  the  circuit  outside  the  gated  region 
will  be  present.  These  effects  tend  to  be  more  pronounced  at  the  band 
edges. 
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D.  Fault  Location 

The  traditional  application  of  the  time  domain  i3  fault  loca¬ 
tion.  This  can  be  done  on  helix  circuits  quite  easily  with  the  ANA. 
Being  able  to  locate  the  general  position  of  a  discontinuity  in  a  TWT 
can  be  a  great  help  in  analyzing  a  possible  failure  due  to  an  opened 
pin-to-helix  connection,  melted  circuit,  etc.  Usually,  a  tube  must  be 
opened  up  to  examine  the  circuit  in  order  to  determine  the  failure 
mechanism.  TDR  can  give  the  user  some  idea  of  what  may  have  happened 
without  opening  the  tube. 

The  best  way  to  implement  this  type  of  use  for  TDR  would  be  to 
plot  the  time  domain  response  of  each  circuit  before  it  is  aged  and 
tested.  Then,  if  a  problem  like  a  melted  helix  is  suspected,  the  time 
domain  response  could  quickly  be  plotted  again,  and  any  significant 
changes  noted.  Although  the  location  of  a  discontinuity  cannot  be 
determined  exactly,  as  was  discussed  in  the  section  on  bead  tracking,  a 
general  id°a  of  its  location  is  obtained.  In  these  types  of  failures, 
the  location  is  not  as  important  as  whether  or  not  they  are  present. 

An  open  pin-to-helix  connection  is  illustrated  in  Fig.  6,  where  it 
is  quite  obvious  what  has  happened,  since  the  response  of  the  circuit 
before  the  open  occurred  is  known.  As  would  be  expected,  there  i3  a 
dramatic  change  in  the  time  domain  response.  An  open  helix  was  not 
actually  tested,  but  was  simulated  by  placing  a  metal  rod  in  the  circuit 
to  create  a  substantial  mismatch.  Figure  7  shows  what  might  be  seen  if 
the  helix  was  melted  somewhere  down  its  length.  From  this  plot,  it 
would  be  easily  seen  that  the  problem  is  not  at  the  input,  as  in  the 
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Fig.  7.  A  melted  (opened)  helix  approximately  one  inch 
past  the  input  connection  in  the  time  domain. 
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IV.  FREQUENCY  DOMAIN  MEASUREMENTS 


A.  General  Theory 


In  this  section,  the  concern  is  to  get  accurate  phase  velocity 


data  utilizing  the  vector  nature  of  the  data  itself  without  having  to 


utilize  the  time  domain  capabilities  of  the  network  analyzer.  This  has 


b°en  done  by  some  researchers  already,  as  reported  in  the  literature, 


however,  the  actual  calculations  involved  are  often  left  to  the  reader. 


It  is  the  purpose  of  the  author  to  include  all  the  necessary  information 


to  allow  the  utilization  of  this  method  to  any  one  with  access  to  a 


vector  ANA  and  a  computer. 


Automation  of  Measurements  and  Calculations 


This  method  is  based  on  a  two-boundary  reflection  model,  which  is' 


discussed  in  Appendix  A.  The  method  of  measurement  is  much  simpler  than 


the  method  used  in  the  time  domain  measurements.  The  first  thing  done 


is  to  measure  the  swept  complex  reflection  of  the  helix  circuit.  At 


this  point,  the  circuits  used  are  actual  circuits  with  loss,  as  this  is 


what  the  model  was  designed  for.  Without  the  bead  in  the  circuit,  it  is 


assumed  that  the  resultant  reflection  is  due  entirely  to  the  input 


discontinuity,  and  that  there  are  no  reflections  further  down  the  tube 


or  that  they  are  hidden  by  the  loss  pattern  on  the  rods.  The  bead  is 


then  introduced,  which  creates  the  second  mismatch  that  is  included  in 


the  model.  The  bead  is  moved  at  regularly  spaced  intervals  along  the 


circuit,  with  swept  data  boing  taken  at  each  position.  These  data  are 


read  from  the  ANA  and  dumped  to  a  floppy  disc  using  a  data  acquisition 


program  on  a  Hewlett  Packard  series  200  PC.  The  data  are  then  read  into 
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an  HP1000  computer  which  has  the  model  program  on  it.  The  whole  process 
could,  of  course,  be  done  completely  on  either  of  the  two  computers. 
Both  programs  are  discussed  in  Appendix  B. 

The  output  of  the  program  at  one  frequency  is  shown  in  Fig.  8. 
The  program  prints  out  the  phase  shift  for  each  position  of  the  bead, 
the  corresponding  phase  constant,  and  the  normalized  phase  velocity. 
Next,  an  average  of  each  of  these  quantities  is  printed  along  with 
normalized  frequency  (ka  cot(ip)),  normalized  phase  velocity 
[  (v/c)cot(i|>)  ],  and  the  high  and  low  values  of  the  phase  shift  and  nor¬ 
malized  phase  velocity.  The  data  were  generally  taken  within  a  uniform 
length  of  the  helix  circuit,  which  is  the  reason  for  averaging  the  data 
over  the  whole  length.  The  data  could,  however,  be  taken  over  a  length 
of  circuit  that  includes  a  pitch  taper,  in  which  case  the  computer  code 
gives  an  indication  of  the  phase  velocity  as  a  function  of  distance 
along  the  circuit.  In  this  case,  the  averaging  of  the  data  at  the  end 
of  the  output  is  meaningless.  The  data  were  taken  in  the  relatively 
loss-free  area  of  the  circuits  tested,  i.e.,  in  the  "clean  length,"  in 
order  to  get  the  best  accuracy  out  of  the  model. 

C.  Comparison  of  Results  with  Known  Phase  Data 


■vXl 
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The  first  circuit  tested  had  a  mean  diameter  of  0.127^  inches  and 
a  pitch  of  0.028  inches,  and  was  designed  for  operation  in  the  2-8  GHz 
band.  This  size  circuit  tended  to  be  the  easiest  to  do  measurements  on, 
as  it  was  much  easier  to  get  the  bead  centered  and  also  to  make  a  bead 
that  was  not  too  large  for  the  circuit.  The  bead  size  used  in  this  case 
was  0.032  inches  wide  with  a  diameter  of  0.0^9  inches.  The  data  were 
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7.40  GHz 


Phi (deg) 

Beta(deg/cm) 

v/c 

107.806 

539.03 

.16529 

104.815 

524.07 

.17001 

103.545 

517.73 

.17209 

106.607 

533.04 

.16715 

109.024 

545.12 

.16345 

109.395 

546.98 

.16289 

110.708 

553.54 

.16096 

105.341 

526.70 

.16916 

111.681 

558.40 

.15956 

107.973 

539.87 

.16504 

103.695 

518.47 

.17184 

107.545 

537.73 

.16569 

107.635 

538.18 

.16555 

Phi  ave.  -  107.367 

Beta  ave.  -  536.83 

v/c  ave.  -  .16605 

kacot(psi)  »  3.2074 

(v/c)cot(psi)  »  .8951 

Phi  varies  by  4.0)1(111.68)  and 

Vrel  varies  by  — 3 . 9X (  .1596)  and 


-3.6*(103.55). 
3- 6%(  .1721). 


Fig.  8.  An  example  of  the  output  from  the 
REFEQ  program  at  one  frequency. 
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Average 

B)°cm) 

Average 

v/c 

346.82 

0.06948 

456.13 

0.06872 

565.99 

0.06815 

682.23 

0.06709 

789.97 

0.06711 

896.06 

0.06721 

1011.0 

0.06671 

1101.5 

0.06782 

1216.7 

0.06736 

1318.8 

0.06769 

1419.3 

0.06816 
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►4.6, 
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►2.7, 

-3.4* 

►3.1, 

-3.6* 

►3.3, 

-4.1* 

►5.4, 

-5.5* 
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constant,  the  normalized  phase  velocity  (v/c),  and  the  amount  of  varia¬ 
tion  in  the  data  for  the  normalized  phase  velocity  at  each  frequency. 
The  bead  was  shifted  at  regular  intervals  of  0.2  cm.  In  Fig.  9,  the 
data  are  normalized  to  the  circuit  diameter  and  pitch  as  normalized  fre¬ 
quency  [ka  cot(i|>)]  and  normalized  phase  velocity  [  (v/c)cot(ij>)  ] ,  where 
cot(tji)  =  it*  diameter/pitch,  and  the  resultant  data  are  plotted  against 
the  data  for  the  same  circuit  obtained  by  the  scalar  bead  pull  method 
described  in  Appendix  C.  As  can  be  seen,  the  agreement  is  quite  good 
over’  the  portion  of  the  band  for  which  comparison  data  were  available. 
Table  7  compares  the  data  at  the  three  points  of  the  vector  ANA  data 
that  fall  close  to  scalar  bead  pull  data  points.  The  data  for  the  two 
methods  are  within  2  percent  of  each  other,  which  is  very  good  agree¬ 
ment  . 


Table  7.  Comparison  of  the  results  obtained  using  the 
ANA  to  calculate  normalized  phase  velocity 
with  that  obtained  using  the  beam  pull  method. 


ka  cot  ( i^) 

HP8510 
( v/c)cot(  \p) 

Bead  Pull 
(v/c)cot(ip) 

Percent 

Diff  erence 

1.0704 

0.9051 

0.8981 

-0.773 

1.3175 

0.8977 

0.8885 

-1.030 

1.5645 

0.8837 

0.8922 

0.962 
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The  next,  circuit  tested  was  the  same  circuit  that  was  tested  in 
the  time  domain  section  for  2-8  GHz.  The  mean  diameter  was  0.151  inches 
and  the  pitch  was  0.088  inches.  The  average  results  of  running  REFEQ  on 
the  data  are  shown  in  Table  8.  The  normalized  phase  data  are  compared 


Table  8.  The  results  of  running  the  program  REFEQ  for  the  second 
circuit  tested.  The  average  phase  constant,  normalized 
phase  velocity,  and  variation  in  the  results  are  shown. 


f (GHz) 

Average 

6 (°cm) 

Average 

v/c 

Variation  from 

Man  in  Data 
for  v/c 

2.0 

137.21 

0.17561 

+3. -5.1$ 

2.6 

179.55 

0.17441 

+3.2,  -3.1$ 

3.2 

222.28 

0.17340 

+3.2,  -3.9$ 

3-8 

266.08 

0.17203 

+3.1,  -4.3$ 

4.4 

311.49 

0.17090 

+12.4,  -9.3$ 

5.0 

352.55 

0.17137 

+10.0,  -8.1$ 

5.6 

398.33 

0.16936 

+4.6,  -3.9$ 

6.2 

443.14 

0.16857 

+4.0,  -3.8$ 

6.8 

487.61 

0.16804 

+4.2,  -3.6$ 

7.4 

536.83 

0.16605 

+3.6,  -3.9$ 

8.0 

579.18 

0.16646 

+4.8,  -3.8$ 

to  the  expected  results  in  Fig.  10.  in  this  case,  the  comparison  data 
were  obtained  using  the  resonance  method,  in  which  case  a  cold  test 
circuit  was  used  with  no  loss  pattern,  as  discussed  in  Appendix  C.  As 
can  be  seen  from  Fig.  10,  the  agreement  looks  very  good.  The  points  are 
compared  numerically  in  Table  9  and  are  seen  to  vary  by  less  than  1 
percent  from  each  other. 


-  32  - 


1.  30 


in 

CL 


O 

U 


/*> 

u 

\ 

> 


KACOT  (PSD 


Fig.  10.  Comparison  of  results  obtained  for  a  second  circuit 
tested  using  the  program  REFEQ  with  those  obtained 
using  the  resonance  method.  Circuit  diameter  and 
pitch  are  0.151"  and  0.088",  respectively. 
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Table  9.  Comparison  of  the  results  obtained  using  the  ANA  to  calculate 
normalized  phase  velocity  with  that  obtained  using  the 
resonance  method  for  the  second  2-8  GHz  circuit  tested 


o 


J 
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IX 


HP8510 
( v/c )  cot  ( 4>) 

Resonance 

(v/c)cot(4>) 

Percent 

Difference 

0.9902 

0.9350 

-0.553 

0.9398 

0.9305 

-0.960 

0.9279 

0.9259 

-0.216 

0.9213 

0.9203 

-0.109 

0.9238 

0.9156 

-0.888 

0.9130 

0.9111 

-0.208 

0.9087 

0.9072 

-0.165 

0.9058 

0.9092 

-0.177 

0.8951 

0.9006 

0.619 

0.8973 

0.8979 

0.067 

ka  cot(i|/) 

0.5635 

0.6935 

0.8235 

0.9535 

1.0836 

1.2136 

1.3936 

1.9737 

1.6037 

1.7337 


The  last  circuit  tested  was  the  IJ  band  circuit  that  was  also 
t°sted  in  the  time  domain  portion  of  this  paper.  The  results  of  calcu¬ 
lating  the  phase  velocity  using  this  method  were  not  much  more  encourag¬ 
ing  than  the  results  found  using  the  gating  method  described  above. 
Figure  11  shows  the  comparison  of  the  data  taken  to  the  expected  data. 
Th°  bead  pull  data  are  what  are  expected  and,  as  can  be  seen  in  the 
figure,  the  data  obtained  here  are  quite  far  off  at  the  high  end. 
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Fig.  11.  Comparisor.  of  results  for  the  IJ  band  circuit 

obtained  using  the  REFEQ  program  with  those  obtained 
using  the  bead  pull  method.  Circuit  mean  diameter 
and  pitch  are  0.0723"  and  0.0*45",  respectively. 
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V.  CONCLUSION 


The  results  of  the  experiments  in  the  time  domain  indicated  that 


precisely  locating  the  position  of  the  bead  in  the  circuit  was  not 


possible.  However,  when  using  the  gating  method  of  determining  a  cir- 


:uit's  phas“  velocity,  good  results  were  obtained  for  the  2-8  GHz  cir¬ 


cuit.  This  indicates  that  although  the  gating  filter  could  not  be 


centered  exactly  on  the  bead's  response,  the  major  portion  of  the  bead's 


response  was  still  included  within  the  filter,  resulting  in  the  expected 


results  for  the  center  of  the  frequency  band.  The  degradation  in  the 


accuracy  of  the  data  obtained  by  this  method  at  the  band  edges  can  be 


attributed  to  this  inability  to  precisely  locate  the  gating  filter  on 


the  bead's  response.  When  attempting  to  utilize  this  method  on  an  IJ 


band  circuit,  no  reasonable  results  were  obtained.  Increasing  the 


resolution  in  the  time  domain  by  increasing  the  frequency  span  of  the 


measurements  did  not  result  in  any  improvement  in  the  data  in  this 


case.  The  tests  were  run  using  both  bandpass  and  the  greater  resolution 


lowpass  modes. 


The  traditional  use  of  TDR,  namely  the  locating  of  faults  as  a 


function  of  distance,  has  proved  to  have  some  use  in  helix  circuits  as  a 


quick  check  on  the  condition  of  a  circuit  in  a  failed  tube  without 


having  to  open  a  vacuum  seal.  Comparing  the  time  domain  response  of  a 


failed  tub^  with  its  original  response  can  show  quite  readily  the  exis¬ 


tence  of  such  failure  mechanisms  as  opened  pin-to-helix  connections  and 


melted  helices. 


The  method  of  obtaining  phase  velocity  by  exploiting  the  complex 
nature  of  the  reflection  data  obtained  on  a  vector  NA  proved  to  be  the 
more  accurate  of  the  two  methods.  Although  both  methods  gave  good 
results  when  used  on  the  2-8  GHz  circuits,  the  frequency  domain  method 
gave  accurate  results  over  the  entire  band  of  interest,  while  the  time 
domain  gating  method  did  not.  The  frequency  domain  method  is  more 
computationally  intensive  than  the  gating  method,  but  would  also  lend 
itself  more  easily  to  total  automation  of  the  measurement  process,  which 
is  a  great  advantage. 

The  failure  of  either  method  to  achieve  any  significant  amount  of 
accuracy  for  the  IJ  band  circuits  tested  was  disappointing.  One  pos¬ 
sible  reason  for  this  is  that  the  size  of  the  circuits  was  much  smaller 
for  the  IJ  band  circuits  than  for  the  other  circuits  tested.  This  made 
it  harder  to  construct  beads  for  these  circuits  and  also  to  get  the  bead 
centered  on  the  axis  of  the  circuit.  Also,  the  match  of  the  IJ  band 
circuit  tested  was  not  as  good  as  the  match  of  the  other  circuits 
tested,  which  causes  the  masking  affects  to  be  more  serious  in  the  time 
domain. 

All  in  all,  it  is  the  opinion  of  the  author  that  the  HP8510  and, 
in  general,  any  time  domain  capable  vector  network  analyzer  can  be  of 
value  in  the  testing  and  design  of  helix  tubes.  The  phase  velocity 
measurement  methods  described  herein  are  quick  and  simple  to  run  once 
they  have  been  implemented  and,  thus,  would  allow  phase  velocity  mea¬ 
surements  to  be  done  on  each  circuit  coming  off  the  line  to  give  some 
idea  how  the  tube  will  behave  in  test.  And,  as  stated  above,  the  ANA  is 


also  of  value  in  determining  the  failure  mechanism  in  a  bad  tube  withou 


having  to  open  the  vacuum  envelope. 

Some  areas  for  further  research  would  be  the  implementation  of  a 
loss  versus  distance  and  impedance  method  using  this  piece  of  equip¬ 
ment.  With  the  addition  of  these  abilities,  the  usefulness  of  the  ANA 
in  the  tube  industry  would  be  enormous. 
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APPENDIX  A 

DERIVATION  OF  PHASE  VELOCITY  MODEL  FOR  FREQUENCY  DOMAIN 

The  model  used  for  the  calculation  of  the  phase  velocity  was  a 
two-boundary  system.  The  circuit  was  assumed  to  have  only  a  single 
discontinuity  at  the  input  to  the  helix  due  to  the  good  match  achieved 
further  down  the  circuit  as  a  result  of  the  loss  pattern.  The  bead, 
when  inserted  into  the  circuit,  is  assumed  to  be  the  only  other  mis¬ 
match.  The  quantity  of  interest  is  then  the  phase  angle  between  these 
two  mismatches.  A  diagram  of  the  model  is  shown  below. 


This  sort  of  system  is  analyzed  in  most  basic  EM  books."  Assuming  that 


r  =  -  r 
2  1 


;he  total  reflection  for  the  system  is 
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Assuming  second-order  terms  to  be  small,  Eq.  A.l  becomes 


r  .  r  +  r  e"J2e 
t  1  3 


Now,  if  one  considers  the  system:  (1)  with  the  second  boundary  removed, 
(2)  with  the  mismatches  an  arbitrary  phase  angle  0  apart,  and  (3)  with 
the  mismatches  an  additional  arbitrary  distance  0^  apart,  the  following 
equations  result: 


r  =  r 

tl  1 


=  r  +  r  e 
til  1  3 


r  =  r  +  r  e 
till  1  ‘3 


-j2(0^02) 


Substituting  Eq.  A. 2  into  Eqs.  A. 3  and  A.  leads  to 


rtII  Ftl  r3  e 


r  -  r  =  r  +  r  e 
tin  ti  i  3 


-J2(e1+e2) 


Eliminating  yields 


rtIII  ~  Ftl  =  ~i2Q' 

FtII  '  rtl  =  6 
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and  from  Eq.  A. 5,  it  can  be  seen  that 
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For  the  exact  case,  one  has 


APPENDIX  B 


COMPUTER  PROGRAMS 


1 .  Use  of  Programs 


a .  Data  Acquisition  Program 


This  program  is  used  to  read  data  from  the  HPS510  and  dump  it  to  a 
floppy  diso  for-  transfer  to  the  HP1000  computer.  The  program  is  in 
basic  and  runs  on  the  series  200  HP  computer.  The  program  prompts  the 
user  for  a  header'  to  write  on  the  first  line  of  each  file  written  on  the 
disc.  It  then  waits  for  the  user  to  press  the  continue  key  when  ready 
to  transfer  data,  i.e.,  when  the  bead  has  been  moved  to  a  new  posi¬ 
tion.  The  program  then  asks  for  the  name  of  the  data  file  to  be 
created.  Each  time  the  sll  data  are  dumped  to  the  disc,  the  data  go 
into  a  new  file.  These  files  must  be  named  using  single  or  double 

letters  in  alphabetical  order,  i.e.,  a,b,c,...,  or  aa,bb,cc .  in 

order  to  use  the  sorting  routine  included  in  the  REFEQ  program.  The 
program  assum°s  a  51  point  trace.  Each  file  on  the  disc  contains  a 
head°r,  the  start  and  stop  frequencies,  and  then  the  trace  data  in 
ascending  order. 


b .  REFEQ  Phase  Velocity  Program 

This  program  is  used  to  analyze  the  data  taken  on  the  ANA  by 
reading  th°  files  created  by  the  data  acquisition  program.  This  code  is 
written  in  Fortran.  It  can  sort  the  data  from  the  floppy  disc  to  create 
the  necessary  input  file,  or  use  an  already  existing  file. 


The  input  file  used  by  REFEQ  is  illustrated  In  Fig.  B.l,  which 


shows  the  necessary  format  of  the  input  file.  The  first  line  is  a 
header  which  can  be  anything,  and  usually  contains  circuit  information, 
date,  etc.  The  second  line  contains  the  mean  circuit  diameter  in  inches 
and  the  value  of  cot(ty).  The  next  lin°  starts  the  first  group  of 
reflection  data  for  a  given  frequency.  Figure  7  shows  a  frequency  of 
2.6  GHz  (first  item  on  line).  The  next  item  on  the  line  is  read  as  the 
variable  N  and  tells  the  program  how  many  times  it  will  calculate  the 
phase  velocity  for  the  given  frequency.  This  number  must  be  one  more 
than  the  desired  number  of  times  the  phase  velocity  is  to  be  calcu¬ 
lated.  In  this  case,  the  19  means  that  the  phase  velocity  will  be 
calculate  18  times.  Also,  there  must  be  N+l  lines  of  reflection  data 
between  this  line  and  the  next  line  containing  an  integer.  Th®  last 
variable  read  on  this  line  is  the  delta  z,  the  amount  that  the  bead  was 
moved  each  time,  in  centimeters.  The  next  line,  line  9,  contains  the 
complex  reflection  of  the  circuit  without  the  bead  present.  The  third 
number  on  this  line  is  an  offset  that  allows  the  user  to  add  a  constant 
multiple  of  180  degrees  to  the  result  of  the  code's  phase  angle  calcula¬ 
tion.  This  is  for  cases  when  the  distance  that  the  bead  is  shifted 
results  in  a  phase  shift  greater  than  180  degrees,  since  the  functions 
used  in  the  code  can  only  give  results  between  ±180  degrees.  Th°  next  N 
lin°s  are  complex  reflection  data  for  the  bead  at  each  subsequent  posi¬ 
tion  followed  by  the  offset.  After  reading  these  values,  the  program 
wants  to  read  N  again,  and  if  it  is  zero,  the  code  is  ready  to  go  to  the 
next  set  of  frequency  data.  If  it  is  nonzero,  the  computer  will  expect 
to  read  N  more  lines  of  reflection  data.  This  is  used  if  the  data  taken 
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6/20/87  0«TS“ OH  611 


.131  3.3907 

“7.6  19  “2 " 
-.41 359E-01 

. 23746E+00 

0 

- ,434 1 9E-01 

. 22393E+00 

0 

'  - . 34878E-01 

.23413E+00 

0  - 

-.4861 5E- 01 

.24855E+00 

0 

-.33783E-01 

•  24741 E+QQ 

0 

- . 29303E-01 

•  2329  0E  +  00 

0  ■ 

-.41 977E-61 

. 22444E+0Q 

0 

- .34481E-01 

.  23339E+00 

0 

-  149037E-01 

“‘.24847E400 

0_  ... 

-  .  33249E- 01 

. 24758E+00 

0 

-  .  29 076E-01 

.23228E+00 

0 

- - . 421 37E-01 

““  .22446E  +  Q0' 

0 - 

- . 54359E-01 

. 23403E+00 

0 

- .  48294E-0T 

.  24839E  +  0  0 

0 

— ",  33524E-01 

"T  24694E  +  00 

0 

- .  29983E- 01 

.  23272E  +  00 

0 

- . 42595E- 01 

. 22472E  +  0  0 

0 

- . 34222E-01 

'  '.23400E  +  00 

0 

-.49156E-01 

.  24873E  +  00 

0 

-. 32448E-01 


. 2474  0E  +  00 


3.2  19  .2 
,21 046E+00 
“  19534E+00 
.21  0  04E+00 
.22467E+00 
“.'21  083E+00 
.  19638E  +  00 

_ .  20870E+00 

■22462E+00 
.21236E+00 
. 1 9576E+  00 
, 2  084  OE+OO 
.  22478E+00 
.21303E+00 


.  1  0998E  +  O  0  0 

.  1 1 162E+00  '0 
.1251 6E+00  0 

.11091 E^O  0  0 

",  95924E-01  '0 
. 1 0899E+00  0 

. 12405E+00  0 

Til 275E+00  "0 
. 9391 7E-01  0 

J^0812E+00  0 

. 12492E+00  0 

. 1 1 273E+00  0 

.93333E-01  0 


TT^9r+«' 
-.  14075E  +  00 
-.  12326E+00 


T16364E  +  00" 
.  14677E-*00 
. 13973E+00 


0“ 

0 

0 


0  0  0  0 


Fig.  B.l.  Example  of  input  data  file  for  REFEQ  program. 
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are  not  evenly  spaced  and  can  be  used  to  exclude  data  that  are  obviously 


bad.  At  the  end  of  the  data  file,  four'  zeros  tell  the  code  there  are  no 
more  data. 

The  program,  when  run,  will  ask  the  user  if  it  is  to  create  an 
input  file  from  data  on  the  disc.  If  so,  the  program  asks  for  the 
number  of  files  to  be  processed  and  whether  they  have  single  letter  or 
double  letter  names.  It  then  prompts  for  the  name  of  the  input  file  to 
be  created,  the  title  line  for  the  file,  the  helix  radius,  and  the  value 
of  co  t  ( vp ) .  If  the  user  does  not  wish  to  create  an  input  file,  but 
wishes  to  us“  an  already  existing  one,  the  code  prompts  for  the  name  of 
the  input  file  to  be  used,  and  whether  the  reflection  data  in  the  file 
are  in  polar  or  cartesian  form.  Next,  the  code  asks  for  the  name  of  the 
output  file  to  be  used  and  the  plot  output  file.  Giving  a  zero  for  the 
plot  file  name  causes  no  plot  file  to  be  used.  Finally,  the  code 
prompts  the  user  to  input  whether  to  use  the  exact  model  or  the  small 
signal  model. 


0 


51 


1 " t; ' I S  I CN  Of-  ccT  05' i  1.3^  R£T4d 

•  this  program  will  retrieve  s  parameter  data  from  the 

'HP  3510  ANA  ANO  WRITE  THAT  OATR  TO  A  DISK  FILE. 

'THE  PROGRAM  ASSUMES  THAT  YOU  ARE  USIN6  A  Si  POINT  CALIBRATION 

A5SI6N  #Hp_ib  TO  7 

ASSIGN  »Nua  TO  716 

AEO°T  9Hp_iO 

CLEAR  #Hp_ib 

REMOTE  @Hp_ip 

REMOTE  9Nua 

DIM  Header*! 601  ,File_neneS!  10]  , BSC  2  3 
DIM  Sll <50 . 1  ) 

OUTPUT  9Nue I "DEBUOFF i * 

OUTPUT  CNwa i "CONT* 

BEEP. 

INPUT  -HEADEP7(TITLE  .  DATE,  ECT..NC1  COMMAS, 30  CHAR  MAX  ,50  DISPLAYED  ON  851 

OUTPUT  9Nwa i ‘ENTOt ‘ i 

OUTPUT  9Nua i ‘ T ITL‘  ‘ ‘  (Header*  r *  ■  i  • 

PRINT  Header* 

BEEF 

OUTPUT  9Nua i ‘ENT0‘ i 
LOCAL  9Nwa 

DISP  ‘PRESS  Continue  WHEN  READY  TO  TRANSFER  DATA" 

PAUSE 

DISP  ‘MEASURIN6  DEUICE* 

REMOTE  9Nwa 

'FIND  the  START  AND  STOP  FREQUENCIES 
OUTPUT  0Nua ( ‘FORMA i  STARi  OUTPACTIi* 

ENTER  CNua i St ar t_freq 

PRINT  Start_freq 

OUTPUT  0N«a i ‘ STOP ( OUTPACT I ; ‘ 

ENTER  @Nuaistop_freq 
PRINT  Stop_freq 

OUTPUT  CNwai  ‘SMi  NUM6I1  F0RM4  (  OUTPDATA‘i 
ENTER  9Nua;S 1  I ( •  ) 

OUTPUT  0Nua i ‘ENTC”  l 
LOCAL  9Nua 
CLEAR  ®Hp_i b 

'STORE  Sll  ON  DISK  WITH  HEADER  AND  FLA6 
'THE  FLA6  ‘T*  SHALL  BE  USED. 

BEEP 

BEEF 

INPUT  ‘FILE  NAME?‘ ,File_nane* 

CREATE  ASCII  FiIe_naneS  ,:5 

ASSI5N  9Path  TO  Fi Ie_name« i FORMAT  ON 

OUTPUT  ®Patni‘T‘ ,Heeder*,Stert_fraq  ,5 1  op_freq ,51 ,SI 1  (  * ) 

ASSIGN  ®Path  TO  • 

PRINT  File_na«eS 
60 TO  40 

INPUT  ‘FILE  NAME7' tFile_naneS 

ASSIGN  ®PatH  TO  F t 1 e_nane* iFORMAT  ON 

ENTER  9PathiB*,Header*,Stert_freq,Stop_freq ,Nnn ,S  I  I ( •  ' 

ASSI6N  CPath  TO  • 

FOR  1-0  TO  50  STEP  5 

PRINT  Start_freq+I»( St op_f req-5 t #r t_freq )/50 , 

FOR  J-0  TO  T 
PRINT  SI  ! ( I  ,J  ). 

NEXT  J 
PRINT 
NEXT  I 

END 
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3.  Phase  Velocity  Code  Listing 


Page  1 


FTN . 


Opt*;  77/LYIS 


Wed  Jun  24,  1937  2:33:08  pn. 

✓DAVIS/PROJ/SAFE.FTN 


1 

2 

3 

4 
3 
6 

7 

8 
9 

1  0 
1  1 
12 
13 
1  4 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 
23 
26 
27 
2B 

29 

30 

31 

32 

33 

34 
33 

36 

37 

38 

39 

40 

41 

42 

43 

44 
43 

46 

47 

48 

49 

30 

31 

32 

33 

34 
33 


%r 


ILC5C  3,3  > 

PROGRAM  REFEQ 


This  program  uses  reflection  date  in  polar  form  (as  per  HP83 1 0  > 
to  solve  for  the  phase  change  between  two  positions  of  the  per¬ 
turbing  object.  The  code  utilizes  a  model  that  assumes  a  single 
mismatch  with  no  perturber  and  a  double  mismatch  when  the  per- 
turber  is  present.  Given  the  complex  reflection  coefficients 
with  no  perturber  and  with  the  perturber  at  two  different  posit¬ 
ions,  the  phase  constant,  the  relative  phase  velocities,  and 
the  normalized  frequency  and  phase  velocity  are  calculated. 


» 

4 


f  req 
mgal , agal 
mga2, aga2 
mga3, aga3 


Input  Variables; 

frequency  at  which  the  group  of  data  points  was 
taken 

magnitude  and  phase  of  the  reflection  with  no  per¬ 
turbing  object 

magnitude  and  phase  of  the  reflection  with  par — 
turber  at  first  position 

magnitude  and  phase  of  the  reflection  with  pci — 
turber  at  second  position 
distance  between  perturber  positions 
allows  user  to  add  an  integer  number  of  pi  to 
phi 

number  of  data  points  in  group.  In  first  line  of 
data  group  for  a  given  frequency,  it  corresponds 
to  the  number  of  computed  angles  minus  one.  Can 
also  be  used  on  a  line  of  its  own  to  restart  for 
non-consecutive  data,  in  which  case  n  more  values 
are  computed,  n  set  to  zero  ends  computation  for 
that  frequency  group. 


a 

4 

4 

4 


delz 
off  set 


4 

4 

4 

4 

4 


Output  Variables; 


phi , phi  a 
vrel, vrela 


beta , betaa 

kacotpsi 
vccotps 1 


phase  difference  and  average  phase  difference  be¬ 
tween  perturber  positions 

relative  phase  velocity  and  average  relative 
phase  velocity 

phase  constant  and  average  phase  constant 
normalized  frequency 

normalized  average  relative  phase  velocity 


REAL  FREQ , HGA1 , MGA2 , MGA3 , AGA 1 . AG A2 , AGA3 , ANG , PH  I , P 1 , DELZ  ,  VREL, 

1  BETA, PHI A, BETAA, VRELA, A, COTPS I .KACOTPSI .VCCOTPS I , RE, IM 

REAL  F, STEP, PCI , PC2 , PC 3 , PC 4 , PH1 1 , PH  1 2, VREL 1 , VREL2 , PH  I MAX, PH  THIN, 
t  VRELMAX, VRELMIN 

IHTEGER  N, OFFSET, NAV.NDF, I, HCHAR 
CHAR ACTER46  0  CASE, HEADER 
CHARACTER* 1 0  OF ILE , IF  I LE , PF ILE 
CHARACTER  MODEL , FLAG , FDAT A , DF ILE42 , FORW43 
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Page  2  REFEQ 


Opt*:  77/LY1S 


Wed  dun  24,  1987  2:38:08  p* 

/DAVIS/PROJ/SAFE .FTN 


36 

57  C 

se  c 

39  C 
60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 
73 

76 

77 

78 

79 

80  1  03 

81  110 
82 

83 

84 
83 
86 

87 

88 

89 

90 

91 

92 

93 

94 
93 

96 

97 

98 

99  113 
1  00 

1  01 
1  02 

103  120 
1  04 
1  03 
1  06 
1  07 
1  08 
1  09 
1  1  0 


COMPLEX  CA1 , GA2, GAS, ARG, DATUMC  2  0, 31  ) 


PI  -  4 . 0*ATANC 1.0) 

ANG  ■  2*P 1/360 . 0 

UR! TE< 1 , *  )  'CONSTRUCT  DATA  FILE? 

READC 1 ,  'C  A  I  )  '  )  FDATA 
IFCFOATA.EQ.  *Y*  )  THEN 

WRITE< 1 ,  •  )  'NUMBER  OF  DATA  FILES  TO  BE  PROCESSED?  ' 
RE ADC  1 >  NDF 


WR I TE< 1  ,  *  )  'NUMBER  OF  CHARACTERS  IN  DATA  FILE  NAMES:  ' 

RE ADC  1.*)  NCHAR 
00  110  I  ■  1 , NDF 

IF( NCHAR ,  EQ .  1  )  DFILE  -  CHAR< 1  +  64  > 

IFC NCHAR . EQ . 2  >  DFILE  •  CHARC 1+64  >//CHARC 1*64  ) 
0PENC8.FILE-DFILE) 

READC8, '<  A1 >' >  FLAG 
READC 8,  'C  A60  )'  >  HEADER 
RE ADC  8 , *  1  FSTART , FSTOP, NPO INTS 
IFCNPOINTS .HE. 31  >  GOTO  300 
00  103  J  •  1 , NPO INTS 
RE ADC  8, *>  RE 
RE ADC  8 ,  *  >  IH 

DATUMC I , J  >  •  CMPLXCRE.IM) 

CL0SEC8) 

UR  I  TEC  1 , * )  'INPUT  FILENAME: 

READC1.2)  IFILE 
0PENC8.FILE-IFILE) 

UR1TEC 1 ,  *  ) 'TITLE  LINE  FOR  DATA  FILEC60  CHAR.  MAX.):' 
READC 1 , 'C  A60  > '  >  CASE 
UR  I  TEC  8,  'C  A60  ) '  )  CASE 

UR! TEC  1 , • >  'ENTER  HELIX  RADIUS  AND  COTCPSI): 

REAOC 1  ,  •  )  A , COTPSI 
UR I  TEC  8 , *  >  A. COTPSI 
URITEC  t , *  >  'ENTER  DELZl  ' 

READC 1.*)  DELZ 
OFFSET  •  0 
F  -  FSTART/1.E9 
STEP  -  3 

DO  120  d  -  1 ,NPOINTS,STEP 
URITEC  8,  •  )  F . NDF-1 , DELZ 
DO  113  I  ■  1 .NDF 

URITEC  8, 116)  REALC DATUMC I, J  )>, AIMAGC DATUMC I , J >  > 
t  , OFFSET 

F  «  F  ♦  $TEP*C C FSTOP -FSTART  VFLOATC NPO INTS” 1  >) 
t  /I  E9 

URITEC8,*)  '0' 

UR  1  TEC  8, *  )  '0  0  0  0' 

FORM  -  'CAR' 

REUINDC  8  ) 


ELSE 


URITEC 1,*)  'INPUT  FILENAME: 
READC 1,2)  IFILE 
OPEHC  8,FILE«IFILE  ) 


-  H8  - 


Page  3  REFEQ 


Opts:  77/LY1S 


1  1  t 
112 

113 

114 
113 
1  16 

117 

118 

119 

120  C 

121  C 

122 

123 

124 

125 

126 

127 

128 

129 

130  10 

131 
1  32 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143  1  1 

144 

145 

146 

147 

148 

149 
130 
lot 
152 
133 
154 
133 

136 

137 

138 

139 
160 
161 
162 

163 

164 
163 


Ued  Jun  24,  1987  2:33:08  pm 

/DAVIS/PRO J/SAFE .FTN 

URITEC 1 , *  )  'INPUT  FILE  IN  <P0L>AR  OF  C  CAR  >TESI AN? 

RE ADC  1  ,  '<A3>'  >  FORM 
END  IF 

URITEC  1 , *  )  'OUTPUT  FILENAME: 

READC 1,2)  OFILE 

URITEC I , *  >  'PLOT  FILENAME l 

READC 1,2)  PFILE 

OPENC7.FILE-OFILE) 

IFCPFILE.NE. '0' >  OPENC 9 , FILE-PF ILE  ) 


URITEC 1 , a  >  'C  E  >XACT  OR  CS>IMPLIFIED  MODEL?  _' 

READC 1 , 'Cal )' )  MODEL 

READC 8,3)  CASE 

URITEC  7,4)  CASE 

IF  C MODEL . EQ . 'E ' )  URITEC7.3) 

IF  < MODEL . NE . 'E * )  URITEC7.6) 

READC  8, *  )  A , COTPS 1 
READC 6, *  >  FREQ,N,DELZ 
URITEC  7 , 13)  FREQ 
NAV  •  0 
PHIA  -  0. 

BETAA  »  0. 

VRELA  -  0, 

PH1 1  -  0. 

PHI 2  -  10000. 

VREL1  -  0. 

VREL2  -  1 . 

READC  8 , *  )  MGA1,AGA1, OFFSET 

IFCFORM.EQ.  'POL' )  GA 1  -  CMPLXC MCA  1 *COSC ACA 1 *ANC >, 

C  MGAl*SINCAGA1*ANG) ) 

IFCFORM.EQ. 'CAR' )  GA1  ■  CMPLXC MCA! , AGA 1  ) 

READC 8,*)  MGA2.AGA2. OFFSET 

IFCFORM.EQ.  'POL' )  GA2  ■  CMPLXC  HGA2*C0S<  AGA2*ANG  >, 

C  MGA2*S INC  AGA2+ANG  )  > 

IFCFORM.EQ. 'CAR' )  GA2  -  CMPLXC MGA2 , AGA2 > 

NAV  ■  HAV  ♦  H  -  1 
DO  20  1-1. N-1 

RE ADC  8 , *  >  MCA3.AGA3, OFFSET 

IFCFORM.EQ. 'POL' )  GAS  -  CMPLXC HGA3*C0SC AGA3-AHG  ), 

C  MCA3*SINCAGA3*ANG)) 

IFCFORM.EQ. 'CAR' )  CA3  -  CMPLXC MGA3, AGA3  ) 

ARC  ■  CGA3-GA1  >/CGA2-GA1 > 

I FC  MODEL . EQ . 'E ' >  ARC  •  ARG-C 1  . 0-GA1*GA2  VC  1 . 0-CA3-GA1  > 
PHI  -  ATAN2C  AIMAGC  ARC ),REALC  ARC  )  VANG 
IFCPHI.LT. 0)  THEN 

PHI  •  -PHI/2.0 

ELSE 

PHI  ■  -PHI/2.0  ♦  180 
EHDIF 

PHI  -  PHI  ♦  180.0  •  OFFSET 
BETA  -  PHI/DELZ 

VREL  •  2E9-PI-FREQ/C BETA-2 .99E1 0-ANG  ) 

PHIMAX  ■  MAXCPHI1 .PHI > 

PHIMIH  •  MINC PHIZ, PHI  ) 


age  4 


REF  EO 


Opt*!  77/LY I S 


Wad  Jun  24.  1987  2:38:08  pm 

/DAVl S/PRO J/SAFE  .  FTH 


VRELMAX  •  MAXC  VREL1  ,  VREL  5 
VRELHIH  -  MIHC VREL2, VREL  ) 

VR I TE< 7.23)  PHI , BETA, VREL 
PHI A  •  PHIA  *  PHI 
BETAA  -  BETAA  ♦  BETA 
VRELA  -  VRELA  ♦  VREL 
PH1 1  -  PHINAX 
PHI2  ■  PHIMIH 
VREL I  -  VRELMAX 
VREL2  •  VRELHIH 
GA2  -  GA3 
RE ADC  8 , *  >  H 
IFCN.NE. 0>  GOTO  11 
PHI  A  -  PHIA/REALC  HAV  ) 

BETAA  -  BETA A/REAL CHAV) 

VRELA  -  VRELA/REALC  HAV ) 

WRI TEC  7,30)  PHIA, BETAA. VRELA 
KACOTPSI  -  2. 0*PI*FREQ*A*C0TPSI/1 1  .8 
VCCOTPS I  -  VRELA*COTPS I 

IFC  PFILE . HE . ' 0 ' )  HRITEC9.3I)  FREQ , PHI  A , BETAA, VRELA , KACOTPS I , 
£  VCCOTPS I 

WRITEC  7,33)  KACOTPSI , VCCOTPSI 
PCI  -  <  PH1 1 -PHIA  >*  1 00/PH1A 
PC2  -  <  PHI2-PHIA  )*1 00/PH I A 
PC 3  -  CVREL1-VRELA)*1 0 0/ VREL A 
PC4  -  CVREL2-VRELA>*10 0/ VREL A 
WRITEC7.36)  PC  1 , PH1 1 , PC2.PHI2 
HR ITEC  7,37)  PC4, VREL2.PC3, VREL1 
REAPC  8 , * )  FREQ.H.PEL2 
IFCCFREQ.HE. 0  )  . AND . < H . HE . 0))  GOTO  10 
STOP 

HR I  TEC  1 ,*>'HP0tHTS  HOT  EQUAL  TO  31!!’ 


FORMAT <  A  1 0  > 

FORMATC  A60  ) 

F0RMATC3X, 'Two  boundary  model  using  simultaneous  aquations  ’ ,/, 

C  7X.A60) 

FORMATC 7X , Exact  modal  used  *****  * ) 

FORMATC 7X, '****•  Small  signal  modal  used  •****•) 

FORMATC ///3X, 'Fraq  -  ' , F3 . 2,  0Hz  ' ,//, 1  OX,  *PhiC dag > ’, 8X, 

C  ’BataC  dag/cm  )  * , 1  OX,  *v/c  ’,/, 1 0X,8C  *- *  >,8X, 12C 1  OX, 

C  3C  ’  ),//> 

FORMATC 1 0X.F8.3,  1  OX.FB.2,1  0X.F7.3) 

F0RMATC/13X, ’Phi  ava .  ■  * ,F8 . 3, /, 13X, 'Bata  ava.  ■  *,F8.2,/, 

£  1 3X ,  'v/c  ava.  •  ' , F8 . 3 ) 

F0RNATC3X.F6.3, 2X, F8 .3,2X,F8.2,2X,F8.3, 2X , F7 . 4 , 2X, F7 . 4  > 

FORMATC 13X,  ’kacotCpsl  )  ■  ' , F7 . 4 , /, 1 3X, 'C v/c  )cotC psi >  *  *,F7.4) 
FORMATC  13X, ’Phi  varies  by  * ,  f3 . 1 ,  #XC  ’ ,  f  6  .2,  *  >  and  *,f3.1,’XC\ 
t  f 6.2, *  >. *  > 

FORMATC 13X.  *Vral  varies  by  ' ,  f3 . 1 , 'XC  ’ , f 6 . 4, *  )  and  ',r3.1,'XC', 
£  f  6 . 4 ,  '  > .  '  > 

FORMATC 2C2X, El  1 .3),2X, 11  ) 

EHO 


mm 


APPENDIX  C 


DESCRIPTION  OF  METHODS  USED  TO  GET  COMPARISON  DATA 


1 .  Resonance  Method 

This  method  of  calculating  the  phase  velocity  of  a  circuit  is 
accomplished  by  building  a  small  length  of  resonant  circuit  into  which 
power  is  loosely  coupled.  The  transmission  characteristics  can  then  be 
examined  on  a  scalar  network  analyzer.  At  each  frequency  for  which 
there  is  an  integral  number  of  half  wavelengths  along  the  total  length 
of  the  circuit,  a  peak  will  be  seen  in  the  transmission.  The  frequen¬ 
cies  at  which  these  peaks  occur  are  recorded,  and  the  phase  velocity  is 
determined  by 


vp  =  2*L*f/N 


where  L  is  the  length  of  the  resonant  circuit,  f  is  the  frequency,  and  N 
is  the  number  of  half  wavelengths  present  at  that  frequency. 


2.  Bead  Pull  Method 

For  this  method  of  phase  velocity  measurement,  a  metal  bead  is 
pulled  through  the  circuit.  The  bead  is  calibrated  so  that  it  moves  a 
known  distance  along  the  circuit.  The  reflection  of  the  circuit  will 
have  a  sinusoidal  type  variation  due  to  the  bead  moving  through  field 
tnaximums  and  minimums.  From  these  data,  the  wavelength  of  the  circuit 
can  be  determined  and  the  phase  velocity  obtained  from 


vp  “  f*wavelength 


where  f  is  the  frequency  at  which  the  data  were  taken.  This  method  can 
be  used  on  either  cold  test  fixtures  or  actual  circuits,  but  whichever 
is  used,  it  is  important  to  have  a  good  match  in  order  to  get  accurate 
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